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Abstract 

The magneto-optical properties of (Ga,Mn)As have been determined within density functional 
theory using the highly precise full-potential linear augmented plane wave (FLAPW) method. A 
detailed investigation of the electronic and magnetic properties in connection to the magneto- 
optic effects is reported. The spectral features of the optical tensor in the 0-10 eV energy range 
are analyzed in terms of the band structure and density of states and the essential role of the 
dipole matrix elements is highlighted by means of Brillouin zone dissection. Using an explicit 
representation of the Kerr angle in terms of real and imaginary parts of the tensor components, a 
careful analysis of the Kerr spectra is also presented. The results of our study can be summarized 
as follows: i) different types of interband transitions do contribute in shaping the conductivity 
tensor; ii) the dipole matrix elements are important in obtaining the correct optical spectra; iii) 
different regions in the irreducible Brillouin zone contribute to the conductivity very differently; 
iv) a minimum in the Re a^x spectra can give rise to a large Kerr rotation angle in the same energy 
region; and v) materials engineering via the magneto- optical Kerr effect is possible provided that 
the electronic structure of the material can be tuned in such a way as to enhance the depth of the 
minima of Re axx- 
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I. INTRODUCTION 



The magneto-optical Kerr effect (MOKE), discovered in 1877, \^ consists in the rotation 
of the polarization plane of linearly polarized light with respect to that of the incident 
light reflected from a magnetic solid surface. But only very recently has this property 
become the subject of intense investigations. The reasons are twofold. First, one can 
exploit this effect to read suitably magnetically stored information using optical means 
in modern hirfi density data storage technology, erasable video and audio disks (magneto- 
optica. di.l.)iQy.fla Second. MOKE is now regarded as a powerful probe in ,nany fields 
of research, such as a microscopy for domain observation, [7| surface magnetism, magnetic 
interlayer coupling in multilayers, plasma resonance effects in thin layers, and structural and 
magnetic anisotropics, [s], ^ Magneto-optical measurements are also a valuable tool in the 
study of the magnetic properties and electronic structure of magnetic materials. Further, it 
is well known that optical reflection measurements can be used to determine the diagonal 
elements of the dielectric tensor. However, the reflection spectra of many intermetallic 
compounds do not show a pronounced flne structure and so the information on the electronic 
structure obtained through reflection spectra is not very detailed. In contrast, the Kerr effect 
shows a finer structure and gives important insights into the properties of transition-metal 



and rare-earth compounds. |lO 
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In the last ten years, transition metal doped semiconductors have attracted considerable 



attention in the field of spintronics, 



ISj where Gai -rMn,. As is 
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jy far the most studied 
Interestingly, the 



material suitable for semiconductor-based spintronic devices. [16 
giant magneto-optical effects observed in some of these systems are especially attractive for 



optical applications. 
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2l| For instance, optical isolators based on magnetic semicon- 



ductors might be ideal components for high speed optical communication systems. 2^ Fur- 
thermore, magneto-optical measurements have been playing an important role in clarifying 
the exchange interaction and the electronic structure of these ferromagnetic semiconductors, 
as they provide very detailed information on the influence of broken time-reversal symmetry 



in itinerant electron quasi-particle states. 
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Despite the increasing role of magnetic semiconductors over the last years, and the exten- 



sive studies available on their electronic and magnetic properties. 
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28|, 



29| 



only very recently have ab-initio calculations of their magneto-optical properties been 
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performed. 



30, ImI, l32| Weng et al. studied the electronic structure and polar magneto-optical 



Kerr effect of transition metal chalcogenides, such as CrSe, CrTe, and VTe in zinc-blende 
and wurtzite structures, by full-potential density-functional calculations 3l|, while in a fol- 
lowing paper Weng et al. analyzed the magneto-optical response of Zni_^.Cr^.Te ordered 



alloys. In Ref. |30|, some of us investigated the magneto-optical properties of Gai.^Mnj-As 
with highly precise first-principles density-functional FLAPW calculations in order to per- 
form an accurate comparison with experiments, focusing on the effect of Mn concentration 
and of the occupied sites. 



The purpose of the present work is to go beyond the analysis presented in Ref. !30|, by 
performing an extensive study focused on the optical conductivity tensor as well as on the 
Kerr spectra. Here we mainly focus on (Ga,Mn)As ferromagnetic semiconductors in the 
high concentration limit by considering a 25% concentration of Mn substituting for Ga. 
Although this can be considered as a rather unrealistic case due to the low solubility of 
Mn in the host GaAs crystal, important aspects of the interplay between the underlying 
electronic structure and the MO Kerr effect (MOKE) can be gained and then extended to 
the low-concentration limit. We therefore studied the optical conductivity in connection 
with the electronic structure. As analysis tools, we considered the density of states (DOS) 
and band structure in order to elucidate the origin of the features of the optical spectra 
in terms of electronic transitions. The origin of these features are further investigated by 
performing an analysis of the electronic transitions throughout the Brillouin zone (k-space 
dissection). In this way, the contributions coming from the different regions of the Brillouin 
zone can be separated out and the role of the dipole matrix elements can be quantitatively 
analyzed. Finally, in order to compare our results with available experiments, we considered 
a more dilute case, namely x=6.25 % . 

After the analysis of the optical properties, we turn our attention to the study of the 
Kerr spectra. Our approach is based on a real representation of the Kerr angle. We show, 
using such a representation, that useful analytical relations can be derived and can be used 
to gain insights into the microscopic quantities that determine the magnitude and frequency 
position of the main features in the Kerr rotation spectrum. 

The work is organized as follows: in the next section, we briefly describe the theoretical 
methods and provide computational details; in Sect. IIIIl we review the electronic and 
magnetic properties of Gao.75Mno.25As; Sect. HV] is devoted to the analysis of the optical con- 
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ductivity tensor in terms of the band structure; in Sect. |V]we perform a full fc-space analysis 
in order to highlight the role of dipole matrix elements in shaping the optical spectra; in 
Sect. IVll we analyze in detail the Kerr spectra; finally, in Sect. IVIIl we draw our conclusions. 



II. THEORETICAL FRAMEWORK AND COMPUTATIONAL DETAILS 



In this work, the Kohn-Sham equations are solved self-consistently, using the full 



potential linearized augmented plane wave (FLAPW) method. 33| We used the local 



spin density approximation (LSDA) for the exchange-correlation functional, with the 
parametrization of Hedin-Lundqvist.[34(] The spin-orbit coupling (SOC) is essential to obtain 



the orbital magnetic moment as well as the magneto-optical effects: 35|, l36|] in the evaluation 
of the optical conductivity, the spin orbit effect is neglected in the self-consistent iterations 
but is included in a second variational step. 

The optical conductivity tensor is calculated according to the Kubo formula in the linear 



response theory: 
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38, 



39| 



n.n' 



(1) 

where a,f3=l, 2, 3; Pa are components of the momentum operator, f^^ is the Fermi dis- 
tribution function ensuring that only transitions from occupied to unoccupied states are 
considered; \kn > is the crystal wave function, corresponding to the Kohn-Sham eigenvalue 
^kn with crystal momentum k; and the delta function warrants total energy conservation. 
The above formula considers infinite lifetime of excited Bloch electronic states. In order 
to take into account finite lifetime effects, we broadened the optical spectra by fixing the 
interband relaxation time for excited states to 0.3 eV. 

Equation [1] contains a double sum over all energy bands, which naturally splits into the 
so-called interband contributions, i.e., n =^ n , and the intraband contributions, n = n , that 
is: 

^.,/^=<r'M+<rM (2) 

For the diagonal tensor components, both terms are important and should be considered 
simultaneously. For metals (or half- metals), the intraband contribution to the diagonal 
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component of a is usually described by the phenomenological expression according to the 
Drude-Sommerfeld model: 40l. |41| 



47r[(^) - luj] 

Within the Drude theory framework, the complex conductivity is fully characterized by 
two parameters: the plasma frequency uop and the relaxation rate 71 = 1/ti. The intraband 
relaxation time, ri, characterizes the scattering of charge carriers, which depends on the 
amount of defects and therefore varies from sample to sample. Here we choose a fixed 
value of 0.7 eV for 71 in all our simulations. It is worth noting that this value does not 
affect our discussion below, and, in any case, variations of its value between 0.2 and 0.7 
eV, which include the usual variability range for 7i,|42)] lead to a negligible effect on the 
optical conductivity for energy values larger than 1 eV. The unscreened plasma frequency 
is obtained integrating over the Fermi surface: 

Up a = ^ I < kns\pi\kns > ^^(ekn - Bp) (4) 

kn,s 

where V is the volume of the primitive cell, ep the Fermi energy, e the electron charge 
and s is the spin. The calculated value of up is 2.42x10^^ Hz (2.75 eV). The intraband 
contribution to the off-diagonal optical conductivity is very small and is usually neglected 



in the case of magnetically ordered materials. 43] 



In the present work, we consider the Kerr effect in the so-called polar geometry, where 
the magnetization vector is oriented perpendicular to the reflective surface and parallel to 
the plane of incidence. In such a case, the Kerr rotation angle 9k{uj) and its elhpticity 77fe(co') 
can be obtained from the conductivity tensor as: 
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'k{ 



'uj)+trik{uj) = , """^^""^ = (5) 



Our study is based on the supercell approach where one of the Ga atoms in an 8-atom 
cell of zinc-blende GaAs is replaced by a Mn atom, thus simulating an ordered alloy with an 
x=25% Mn concentration. The lattice parameter is chosen equal to the experimental lattice 
parameter of GaAs (5.65 A). For larger unit-cells, i.e. x = 6.25 %, we refer to the work of 
Picozzi et al.,^^ where the pertinent computational and structural details are given. 



As far as the other technical details are concerned, we used a wave-vector cutoff of the 
basis set equal to Kmax=3-5 a.u. and an angular momentum expansion up to Imax = 8 for 
both the potential and charge density. The muffin-tin radius, Rmt, for Mn, Ga and As were 
chosen equal to 2.1, 2.3, and 2.1 a.u., respectively. The relaxed internal atomic positions 
were obtained by total-energy and atomic-force minimization using the Hellmann-Feynman 
theorem 4J| with residual forces below 0.015 eV/A. The sampling of the irreducible wedge 



of the Brillouin zone (IBZ) was performed using the special fc-point method. 45|] In order 
to speed up the convergence of the sampling, each eigenvalue is smeared with a Gaussian 
function of width 0.03 eV . The ground-state electronic structure is calculated using a 
4x4x4 cubic mesh. We checked that these computational parameters are accurate enough 
to obtain total energies and magnetic moments within 10-15 meV/Mn and 0.01 fiB (keeping 
the muffin-tin radii fixed), respectively. On the other hand, in order to accurately compute 
the optical conductivity, the k integration must be carefully taken care of. We investigated 
its accuracy by varying the number of k points up to a (10,10,10) shell (35 fc-points in the 
IBZ): we get the Kerr rotation angle and ellipticity within a few hundredths of a degree 
already using a (4,4,4) shell. 



III. ELECTRONIC AND MAGNETIC PROPERTIES 



Before discussing the magneto-optical (MO) properties, we consider in some detail the 
electronic structure and the magnetic properties of Gan 75 Mnn -^.^ As . Although shown and 
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52 



53| we show them in 



discussed in several previous publications. 
Fig. [1] and Fig. [21 as it is fundamental to have in mind a clear picture of the electronic and 
magnetic properties of the compound considered. This is needed in order to consistently 
follow the discussion of the optical conductivity tensor (see below), which is the main issue 
of the present work. In particular, in Fig. [1] we show the spin polarized total density of 
states (DOS) and the site angular momentum-projected DOS (PDOS) for Mn and nearest- 
and next-nearest-neighboring atoms, that is As and Ga atoms, respectively. In Fig. [2l we 
show the spin-polarized band structure (not including spin-orbit coupling), along the F-X- 
M-F-R-M symmetry directions, for the majority (left) and minority (right) bands. Numbers 
in round brackets label groups of bands (eventually degenerate) at the F point. Vertical 



7 



dashed lines denote the electronic transitions giving rise to spectral features in Re Oxx (see 
next Section), and in the following, we discuss both Figs. [1] and [2l The total DOS (see 
FiglH topmost panel) shows the marly half-metallic (n-HM) nature of the electronic ground 
state: the majority bands cross the Fermi level (set to zero) whereas a net gap opens in the 
minority bands with the Fermi level cutting the minority conduction band near F. 

The total DOS curve is composed of several characteristic structures whose origin can 
be easily recognized by looking at the partial (P)DOS curves. In the following, we only 
consider the upper part of the valence band (above the ionicity gap, that is above —7.5 eV) 
and the lower part of the conduction band (up to ~ 5 eV) which define the energy window 
relevant for the optical transitions we are interested in. The peak centered at —6 eV in 
the DOS curve (for both spin components) is dominated by the energetically low-lying Ga-s 
states, and corresponds to the lowest lying bands in Fig. [2], [(1) and (la)], which are rather 
(surprisingly) well dispersed along F-X-M-F, and almost fiat in the F-R-M direction. As the 
energy increases, there is a peak at —5.0 (—4.6) eV in the majority (minority) component, 
which mainly derives from s-Mn; it corresponds to the next low-lying band shown in Fig. [21 
[(2) and (2a)], which is degenerate with the lowest band [(1) and (la)] at R. 

A broader structure in the total DOS extends from about —4.0 eV up to Ep (majority 
component) and up to about —1.0 eV (minority component). Let us consider the majority 
component. The sharp peak at —3.0 eV originates from Mn-ci orbitals and corresponds to 
the group of bands in Fig. [2] [(3)]; these bands are degenerate at F and are split into three 
bands along F-M and into two bands along F-R. Inspection of the eigenvalue decomposition 
shows that these states have ^2 symmetry at F and thus they hybridize with anion p states 
(see the corresponding PDOS curve in Fig. [1]). In general, the higher energy bands [groups 
(4) and (5) in Fig. [2] mostly derive from hybridization between Mn-c? and As-p states, as well 
as the bands [(6) and (7)] around Ep. The Mn-s states contribute little to the defect bands 
since they obviously lie lower in energy than the Mn-d and -p states. One also recognizes 
the intra-atomic hybridization (among states of the same atom) occurring for Mn-p and -d 
states: the corresponding DOS curves overlap in a wide energy range. Higher bands [(8)] 
are mostly derived from cation and anion states. For the minority component, the valence 
bands, [(3a), (4a) and (5a) in Fig. [2], are mostly dominated by As-p states. The Mn-c? 
states are pushed upward in energy above Ep, due to exchange splitting: these states are 
almost unoccupied and remain fairly narrow in energy above the gap. The conduction bands 
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[(6a), (7a) and (8a)] are composed of Mn-d and -p states just above Ep, while higher bands 
[(9a)] mainly by Ga-s and As-p states. The same Mn-p-d intra-atomic mixing is found in 
the minority component above Ep (the Mn-p and -d DOS are localized in the same energy 
range) . 



IV. OPTICAL CONDUCTIVITY 
A. General considerations 



In Fig. [31 we show the (a) real (left) and (b) imaginary (right) parts of the diagonal 
(top) and off-diagonal (bottom) components of the conductivity, in the range 0-10 eV. For 
simplicity, we did not include the Drude term in the conductivity curves shown. In the real 
part of the diagonal component (Re axx), we highlight some features with vertical lines, 
which will be discussed in closer detail later. 

The Re axx curve has a quite broad peak at around 3.8 eV whereas the Im a^x curve shows 
a deep minimum in the low energy range (around 2.3 eV) and a shoulder above 6.8 eV. The 
Re axy curve shows small amplitude oscillations in the energy range 0—4 eV; a deep minimum 
is found at ~ 5 eV, and a maximum follows at ~ 6.8 eV followed again by oscillations with 
small amplitude. The imaginary component, Im axy, has similar characteristics in the low 
and high energy range (small amplitude oscillations), but a maximum develops at ~ 4.6 eV 
and a minimum follows at ~ 6.4 eV. 

The absorptive part of the off-diagonal optical conductivity, Im axy, has a direct physical 
interpretation. It is proportional to the difference in absorption rate of the left and right 
circularly polarized (LCP and RCP) light, and its sign is directly related to the spin polar- 
ization of the states responsible for the interband transitions producing the structures in the 
spectrum. 5^ The peak in Im axy around 4.6 eV suggests that interband transitions related 
to LCP light should be stronger in this energy region. On the other hand, around 6.4 eV, Im 
axy shows a deep minimum value, indicating the dominance of interband transitions related 
to RCP. Finally, Im axy is negligible when the absorption coefficient for RCP is equal to 
that of LCP - i.e., in the lowest and highest parts of the spectrum. 
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B. Analysis in terms of band structure 

Since Re axx is directly related to the density of states and transition probabilities, [5^ 
we investigate the calculated energy band structure in order to analyze the details of the 
spectra. It has been known for quite some time that both the spin-orbit (SO) interaction and 



the exchange splitting are needed to produce a non-zero Kerr effect. j35|, l55|, |56|] Therefore, 
a proper analysis should be performed considering the band structure with spin-orbit (SO) 
included. However, we remark that the calculated band structure with SO (not shown here) 
is quite similar to that formed by the superposition of the majority and minority spin bands 
(without SO). This corresponds to a spin-orbit coupling interaction that is small compared 
to the exchange splitting, so we can qualitatively consider the former as a perturbation 



to the latter: |57| this allows us to restrict our discussion to the usual spin-polarized band 
structure, thus neglecting the mixing of majority and minority spin components. In this 
way, a decomposition with respect to each separate spin channel (spin up, down) can be 
done; this decomposition is almost uniquely determined if the spin-orbit coupling is weak. 
Inspection of the eigenvalue decomposition at V shows that the exchange-splitting of As-p 
and Mn-figg states are ~ 0.03 and 2.5 eV, respectively, whereas the spin-orbit splitting is ~ 
0.02 and 0.002 eV, respectively. In this case, the main effect of the SO interaction is just to 



58|, 



59 



60| 



remove some accidental and systematic degeneracies 

In Fig. [21 we can recognize three different contributions coming from the band structure 
which build up the optical spectra: (a) close parallel bands crossing Ep] (b) bands which are 
degenerate (or almost degenerate) at one point and then separate out, moving below and/or 
above Ep; and (c) occupied /unoccupied bands which are almost parallel in a significant 



part of the Brillouin zone. 6l| The first type does contribute to the conductivity in the low 
energy range and gives rise to bumps of interband transitions; the second type gives rise to 
an almost constant amount of interband transitions with a distinct onset energy; the third 
type gives rise to peaks at higher energies. In the following, we refer to hands although 
a more appropriate discussion should be done in terms of constant eigenvalue surfaces in 
reciprocal space. j62] However, the purpose of this section is to analyze on a qualitative level 
the origin of the different structures in Re axx- For this purpose, we highlight the relevant 
electronic transitions (see Fig. [2]) that give rise to the spectral features in the Re axx curve 
at 0.4, 1.1, 2.8 and 3.8 eV (see Fig. [3]), respectively. 
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Re a XX is non-negligible already at rather low energies (as low as 0.5 eV, see Fig. [3]). As 
can be seen from Fig. [21 interband transitions are available at this energy due to bands which 
are degenerate at some k points at energies slightly below (or above) Ep and then separate 
out, moving above Ep [type (b) transitions]. For example, in the minority component, the 
main contribution in this energy range comes from the two bands [(6a)] below Ep^ which 
are degenerate at F and split along F-X and F-M. Incidentally, we note that the same bands 
are responsible for the n-HM character of the (Ga,Mn)As compound. At very low energy, 
the bands along F-M do contribute to the conductivity, whereas those along F-X do so at 
higher energies. 

Consider in detail the minority bands along F-X with energies just below Ep [bands (6a) 
at F, see Fig. [2]. A careful inspection of the eigenvalue decomposition into atomic-site- 
projected wave functions shows that these bands have a prevalent d character but with a 
non-negligible contribution from anion "p as well as s cation states. Moving away from the 
zone center, the d as well as p contributions increase in the higher energies of these bands, 
while the s contribution decreases; the opposite occurs for the lower band. The d — p and 
s —p mixing is crucial for the onset of the dipole transitions among these states at very low 
excitation energies which otherwise would be forbidden by selection rules. Moving away from 
F, their energy separation increases, as is typical for type (b); as the lower band crosses Ep^ 
the interband transitions stop contributing at an excitation energy around 0.40 eV. This 
corresponds to the first fiat region in Re axx starting around 0.40 eV. Other electronic 
transitions at 0.40 eV are shown in Fig. [2] with lines near the M point in the spin-up bands. 

Upon further increase of the excitation energy, additional type (b) structures start con- 
tributing, like the majority group of bands [(7)] above Ep and degenerate at R which split 
along R-F and R-M. In particular, along R-F, one band remains above Ep, while the other 
goes below. This causes an increase of Re axx-, at energies lower than 1.0 eV. In this energy 
range, type (a) bands also contribute, like the almost parallel majority bands which cross 
Ep along F-M, and those crossing Ep along F-X [all of them deriving from bands (6) and 
(7)]. 

As the photon energy reaches ~ 1.1 eV, several transitions are activated, accounting for 
the bump in the conductivity. The transitions at 1.1 eV are shown with blue (dotted) lines 
in Fig. [21 In the majority (minority) component, they are mainly localized near the R (F) 
point [mostly type (c) transitions]. 
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At ~ 2.8 eV, we have a bump in the Re axx (see Fig. [3]), First, let's consider the minority 
component. The almost dispersionless occupied band around —1.0 eV [(5a)] is coupled to 
the group of bands around 1.5 eV [(8a)]. These are type (c) transitions and they give rise 
to the bump in Re a^x around 2.8 eV. The transitions which match the energy of 2.8 eV 
are shown by red (dashed) lines in Fig. [21 In the majority component there is also a 
contribution coming from the group of bands (5) and (7). As the photon energy increases, 
Re axx reaches its maximum value, at around 3.8 eV. There are several contributions to 
this peak: they involve degenerate minority bands at high symmetry points, like the groups 
of bands at F centered at —2.5 eV [(4a)] and 1.5 eV [(8a)] and at R with energies ~ 
—2.0 eV and 2.0 eV, respectively. In the majority components, more type (c) transitions 
contribute to the maximum of Re axx- They mainly involve groups of degenerate bands at 
F (occupied/unoccupied) which split along F-M: they are not flat along this symmetry line 
but are almost parallel near F. For instance, the highest of this group of bands [(5)] split 
along F-M couples with the lowest of the group of bands [(8)]. 

As the photon energy increases, the lowest majority and minority energy bands between 
—5 and —7 eV below Ep come into play. The dispersion of the occupied as well as the unoc- 
cupied bands involved in the transitions increases, and this correlates with the conductivity 
decrease. 

It can be useful to consider the role of the electrons with spin up and spin down in 
shaping the spectra of Re axx- Again, this can be done on a qualitative level only. In 
fact, the dipole matrix elements enter quadratically in the calculation of the conductivity 
(see Eq. [J]) and interference effects can be lost by separating the spin-up or spin-down 
contributions to the conductivity. Indeed, from Eq. [H we have, for a=(3=x and neglecting 
the k index: 



< n\px\'n >< n \px\n >= \ < n\px\n > |^ (6) 
where n includes the spin indices. Expressing the spin indices, and neglecting spin-flip 



transitions 



63 



6J| we have: 



< n\px\n > 



P 



I n |2 + ipU ,|2 + 2/?e(p^T , 

' x.nn ' ' x,nn ' ^ x,nn x,nn 



(7) 
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In Fig. m we show the two spin contributions to Re axx-, i-^- neglecting the last term 
in Eq. [3, in the range 0-10 eV. Incidentally, we note that the sum of the two separate 
contributions (which does not include the interference effects coming from the last term in 
Eq. [7]), not shown in Fig. HJ differs in some fine details from the curve shown in Fig. [3] 
(which naturally includes the interference effects). This confirms a posteriori that the loss 
of interference effects due to neglecting the last term in Eq. [7] introduces a negligible error 
and allows us to treat the two spin contributions to the conductivity separately. The two 
curves are qualitatively similar, and they clearly show the origin of the bump at 2.8 eV and 
the maximum at 3.8 eV: the former is mainly due to transitions involving minority while 
the latter the majority spin bands. From Fig. [2] one could conclude that both the features 
at 2.8 and 3.8 eV arise from minority spin band transitions: this is not exactly true due to 
the role of the dipole matrix elements. 

In the same Fig. IH we also show the joint density of states (JDOS), which can be derived 
from Eq. [T]by replacing the matrix elements with a constant factor set to 1 in our case.js^ 
This quantity clearly highlights the role of the dipole matrix elements: as a matter of fact, 
the JDOS does not show any relevant feature in the spectra stressing that the spectral 
features observed occur primarily due to the very large matrix elements rather than to a 
large joint density of states effect. On the same basis, the large shoulder above 6 eV in 
the JDOS does not produce any peak in the conductivity spectra due to the small matrix 
elements. 



C. Comparison with experiments 



Finally, we compare our calculations with some experimental results. For this purpose, 
we consider a lower concentration x = 6.25% (1 Mn in a 32-atom cell) that is closer to 



available experiments. 65|, l66|, l67|| In Fig. [5] (a) we show the calculated real component of the 



diagonal optical conductivity. Re function of the photon frequency, compared with 

spectra obtained from experiments and model calculations [68| for Gao.95Mno.05As (at hole 
concentration p ~ 0.8 nm~^), in the [0.01-2] eV energy range. In this interval, our calculated 
conductivity is of the order of 1-2-10'^ (ohm-crn)"^, i.e., remarkably larger (almost by an order 
of magnitude) than the experimental values 651] but of the same order of magnitude as the 
results of model calculations. 681] The disagreement with experiments may be fully explained 
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in terms of the suggested values of carrier concentration in the experimental samples: p ~ 
0.35 nm^^ for a Mn concentration on the order of ~6%, with a degree of compensation on 
the order of 70-80%. 

In our systems, a pretty naive evaluation of the hole concentration can be obtained by 
supposing that each Mn - with a nominal valence of +2 - substituting a Ga atom gives rise 
to a hole; therefore, 1 Mn atom in a 32 atom cell produces a hole concentration of 1.38 nm~^. 



It has been shown 68|] that the hole concentration has a relevant influence on the optical 
conductivity: for example, at very low energies (<0.01 eV) cr^a; ranged 68j from ~ 1.5x10^ 
to ~ 3.5x10^ (ohm ■cm) for hole concentrations varying from p ~ 0.03 to 0.80 nm . These 
values are very consistent with the range in which our conductivity falls. Looking at the 
trend as a function of energy, our a^x spectrum is generally featureless, with the exception 
of the minimum at 800-1000 meV - in remarkably close agreement with experimental as 
well as model calculation results. 68| In the inset, we show the effect of including the Drude 



contribution to the spectra in an extended energy range: as expected, only in the low 
energy range (<~ 1.5 eV) does the Drude term significantly contribute to the total spectra. 
In particular, the above mentioned minimum is present in both total and interband-only 
terms, but it is much more enhanced upon inclusion of the Drude term, therefore improving 
the agreement with experiment. 



The real and imaginary parts of the dielectric constant, 



69|, are shown in Fig. 



151(b) 



and (c), respectively, and compared with experimental ellipsometry measurements. 66|] The 



agreement is rather good, as far as the overall features are concerned. However, there is a 
sizeable discrepancy in the energy position of peaks and valleys. This has probably to be 
ascribed to the neglect of self-energy corrections or to the incorrect treatment of correlation 
effects, resulting in an only partial description of the underlying electronic structure. In 
order to improve the agreement with experiment, we used the so-called A- fitting procedure, 
suggested by Rhee et al. (see Ref. [70| for details), which is an oversimplified approach 
to include a rescaling in the excitation energy spectrum, avoiding the complicated task of 
properly evaluating self-energy effects. Using a value of A ~ —0.1 eV, both the real and 
imaginary spectra were reasonably reproduced. 
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V. k SPACE DISSECTION 



In this section, we study the origin of the spectral features of Re axx- The quahtative 
analysis discussed in the previous section is based on the band structure which highlights 
only transitions involving electronic states at symmetry points or lines in the irreducible BZ 
(IBZ). Here we examine the spectrum in further detail, taking into account the contributions 
coming from the full IBZ. Thus, we sample the BZ with an 8x8x8 cubic shell and evaluate 
the matrix elements in Eq. [1] in a small energy window (5e=0.2 eV) centered around the 
spectral feature of interest. For each k-point, the contribution to the conductivity {c^) 
was evaluated calculating the dipole matrix elements between the initial and final states 
contributing within the selected energy window. Finally, the have been normalized to 
unity, . At each k-point, the is graphically represented by an arrow whose length is 
proportional to its magnitude. The results are shown in Fig. [61 

From this plot, one can easily visualize the region of k-space contributing to the features 
in the conductivity spectra. We focus on the maximum of Re axx-, at 3.8 eV. Figure [6] shows 
the origin of all interband transitions in the range 3.6-4.0 eV, contributing to the maximum 
of Re axx- First, we note that different regions in the IBZ contribute to the conductivity 
very differently: the largest contributions come from the region around the R point, and 
the strength of the contributions decreases upon moving towards M and X. The k-points 
where large interband transitions occur belong to a line parallel to M-R, suggesting that 
these are transitions between bands or band pairs which are very flat throughout much of 
this symmetry line or transitions with very large dipole matrix elements. The contributions 
from planes parallel to R-X-M rapidly decrease in amplitude and become negligible for k- 
points close to F. From this, it is clear that the most active region in the IBZ in shaping 
the maximum of Re axx is the region parallel to the R-X-M plane, near the edge of IBZ. 

VI. MAGNETO-OPTICAL PROPERTIES 

We now turn our attention to the Kerr spectra. Equation [5] suggests that the Kerr an- 
gle could be enhanced by a large magneto-optical component {axy) and a small optical 
component {axx)- To investigate the Kerr spectra on a qualitative level, one can con- 
sider the separate contributions of the numerator axy{ijj) and the denominator D{ijj) = 
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<yxx{'^)\h- + cTxxi^)- The corresponding features in the imaginary part of the spectra are 
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then correlated to those observed in the Kerr spectra. This analysis has been done in Ref. 
by some of us, to which we refer the reader for further details. Incidentally, we note that our 
spectrum (presented in Fig. [71 see below) is quite similar to those presented in Ref. [sol for 
X = 6.25% and x = 12.5%. We can conclude that the main features of the Kerr spectra of 
(Ga,Mn)As do not depend to a great extent on the Mn concentration, at least for the cases 
considered: this further validates our choice to focus on the high concentration limit in this 
study. 

As pointed out by Schoenes and Reim, {tiJ , the form given by Eq. [5] of the Kerr rotation 
is not the most suitable to discuss due to the complex denominator. Indeed, it is possibile 
to express the real part of the left hand side of Eq. [5] (that is the Kerr angle) in terms of the 
real part of the complex fraction appearing in the right hand side of Eq. [5l In this way, one 
ends up with a real fraction where the denominator ("optical component") depends only on 
Re axx and Im axx and the numerator (" magneto-optical component" ) depends on both the 
real and imaginary parts of a^x o^nd the Appendix): 

with A = RaxyRcTxx — Icrxylo'xx and B = laxyRcTxx — R<^xyl<^xx- (Here, R and I refer to 
Real and Imaginary part, respectively). For the definitions of p, g, Dq, and Di we refer 
to the Appendix. We further define Ni = — A{ijj)p{uj) , N2 = —B{uj)q{uj), N=Ni+N2 and 
D=DoDi so that 6k = N/D. We will show that this decomposition can be useful to study 
the correlations between features of the Kerr spectra and electronic structure. 

In Fig. [TJ we show: ellipticity (sk), and Kerr angle {9k, in decimal degree), in panel (a); 
Ni, N2 and N in panel (b); Dq, Dq ^, and Re axx (the spectrum has been rescaled in such a 
way that Re axx and Dq have the same maximum value) in panel (c); Di, and D^^^ in panel 
(d); and finally D, Re axx (the spectrum has again been rescaled in such a way that Re axx 
and D have the same maximum value), and D~^, in panel (e). The notation means 
l/D... (reciprocal function). All the quantities, except the Kerr angle, are in arbitrary units. 

Let us focus on the Kerr ellipticity and rotation [panel (a)]. Clearly, they are related: 
when the Kerr ellipticity crosses the zero line, a peak appears in the Kerr rotation spectra due 



to the Kramers-Kroning relations. 3l|| It is interesting to note that the theoretical ellipticity 
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curve crosses/touches the zero axis a number of times, hence suggesting that the incident 
hnearly polarized hght, at these frequencies, would stay as linearly polarized only. The 
Kerr rotation is characterized by several spikes as well as sign reversals in the whole energy 
window. In particular, there is a first magneto-optical resonance (0.56°) at /Lj ~ 0.2 eV 
and other main peaks are located at ~ 1, 2 and 6 eV. The arrows in the Kerr spectra mark 
the main features (positive and negative). From Fig. [7] [panels (a),(b)], we see that the 
shape of the Kerr spectra is mainly determined by Ni in the low energy range (0-3 eV), 
by both Ni and N2 in the medium energy range (3-7 eV), and by N2 in the high energy 
range (7-10 eV). Very interestingly, the overall spectral trend of the Kerr angle, such as 
the peaks and sign reversal positions, are very close to that of the numerator N{uj) [panel 

(b) ], and only the relative heights of the spikes are different. As expected, the zero of N{uj) 
fixes the zero in the Kerr spectra. We can draw the first conclusion: N{uj) determines the 
overall trend of the Kerr spectra, and, in particular, the presence of maxima and minima. 
We also stress that the analytic form of the numerator of the Kerr angle (see Eq. [8] and the 
Appendix) entangles in a rather complicated way all the real and imaginary components of 
the conductivity tensor: there is no simple guideline to link the trend of a part of the Kerr 
spectra to the spectral trend of a specific (real or imaginary) tensor component in the same 
energy range. 

We now focus on the denominator D~^. First we write ^x=N(u;)D~^(ci;), so that 
plays the role of an enhancement factor. It is evident that both Dq and Di are positive 
definite, that is they never cross the zero line: the enhancement factor does not have any 
"resonant" behavior. In Fig. [7| panel (c), we can see that, remarkably, the spectral trend of 
Dq is very similar to Re a^x- We might therefore infer that Dq ^ has large values whenever 
Re axx has its smallest ones. In our case, this happens when the photon energy approaches 
~ eV. On the other hand, Di [panel (d)] does not show any remarkable features, apart 
from the maximum at very low energy. Considering that the reciprocal function Dq ^ [panel 

(c) ] exhibits a maximum at ~ eV and strongly decreases with increasing energy, D^^ has 
a minimum at ~ eV and monotonically increases with energy. D^^ is also featureless, in 
almost the whole energy window. From panel (e), we see that D~^ has a large peak around ~ 
0.5 eV, then decreases very fast to an almost constant value, and slowly increases at energies 
larger than ~ 5 eV. Clearly, the shoulder of D~^ at low energy comes from the behavior of 
Dq ^ and in the same energy range, whereas the increase of D~^ at high energy is mainly 
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due to the behavior of D^^ in the corresponding energy range. 

These results clearly explain the origin of the Kerr spike (0.56°) at ~ 0.2 eV. The nu- 
merator determines the existence and location of the spikes: whether they give rise to a 
strong Kerr angle follows entirely from the enhancement factor, D^^, which, in this case, 
magnifies the peaks at low energy and suppresses those at higher energy. The importance 
of is highlighted by again considering N(ct;): we would expect large Kerr angle values 
at higher energy (around 6 eV). However, the corresponding peaks are suppressed by the 
D^^ factor. Hence, we are ready to conclude that the large Kerr peak at low energy has an 
"optical" origin; the features at high energy have an MO origin (they follow the trend of the 
numerator) but they do not result in large peaks due to the damped behavior of in this 
same energy range. This confirms in a more transparent and direct way what was already 



observed in Ref. 



301 ] More important, the similar trend of Dq{uj) (D) and Ra^x suggests 



that the presence of minima in the Raxx spectra may give rise to strong Kerr angles in the 
same energy range. On the other hand, RcTxx is directly linked to the optical transitions, 
that is to details of the electronic structure of the materials: this could have far reaching 
consequences in magneto- optical Kerr effect engineering, provided that the electronic struc- 
ture of the material can be tuned correctly in such a way as to enhance the features of the 
Re a XX minima. 



VII. CONCLUSION 



We presented results of first-principles calculations of the magneto-optical properties of 
(Ga,Mn)As within density functional theory aimed at investigating in great detail the role 
of the electronic and magnetic properties in determining the magnetooptical behaviour of 
the material. The spectral features of the optical tensor in the 0-10 eV energy range were 
analyzed in terms of the band structure and density of states and the role of the dipole 
matrix elements was highlighted in terms of Brillouin zone dissection. 

We found that different types of interband transitions contribute in shaping the conduc- 
tivity tensor. The dipole matrix elements play a key role greatly affecting the optical spectra 
in the low as well as high energy ranges. Brillouin zone dissection reveals that different re- 
gions in the Brillouin zone (not necessarily symmetry points or lines) can contribute very 
differently in shaping the conductivity tensor; we find that the most active region in the IBZ 
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is around the R point. 

Moreover, since the Kerr rotation spectrum is a resuh of a complex entanglement of real 
and imaginary contrbutions, we presented a possible way to analyze the calculated Kerr 
spectra in terms of a real representation of the Kerr angle. To the best of our knowledge, 
this has not been used as an analysis tool in the past, although it is implicit in Eq. O Using 
this representation, we have clearly elucidated the origin of the features of the Kerr spectra 
and have highlighted the role of the minima of Re axx which can possibly correlate with a 
large Kerr angle (in our case, at low energies). Nevertheless, the results given above support 
the conclusion that at present it is difficult to give simple rules for the occurrence of large 
spectral features at specific laser-light frequencies. Indeed, both spin-orbit and exchange 
splitting effects are entangled in N{uj), which determines the overall behavior of the Kerr 
spectra and eventually the strength of the Kerr features. No simple rules can be given to 
disentangle their effects on the absolute magnitude of the Kerr angle. On the other hand, 
we found that part of the Kerr spectrum (in our case, the low energy range) is very sensitive 
to the shape of Re axx, which is very closely related to details of the band structure of the 
materials. 

While it is obvious that without a precise knowledge of the material band structure it is 
not possible to make a priori predictions of the Kerr angle magnitude, it is also true that 
a proper tuning of the band structure would allow one to consequently tune the magnitude 
and location of spectral peaks through the Dq term, whose spectral behavior is very close 
to Re axx- However, further studies are needed in order to better elucidate how this proper 
tuning of the underlying electronic structure can be achieved (i.e. by making use of alloying, 
pressure, strain field) in order to enhance the minima of Re axx- Therefore, we expect that 
computational materials design can substantially contribute to magneto- optical Kerr effect 
engineering in the future. 
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APPENDIX 

Here we give an exphcit representation of the complex Kerr angle {6k + ir/k, Eq. in the 
form a+if3, where a and (3 are real functions: the Kerr angle {6^) corresponds to a. To 
this end, we need to calculate the square root of the complex number in the denominator of 



Eq. [51 It can be helpful to make use of the following formula 72] : given \/a + ib = p + iq we 
have 

sgn{b) 



p = -j=\/ Va^+lP + a q= -^^-^ \J ^ a? + - a (9) 



where sgn{b) is the sign of b. In our case, 



1 H c^xx = \/ (1 ^(^xx> + i — R(^xx 

UJ \ iO 00 

where R and I stand for real and imaginary part and a = 1 — ^la^x, b = ^Ra^x- By 
definition, b is positive. 

In order to get rid of the complex number in the denominator of Eq. [5l we multiply 
numerator and denominator by cr^xiP ~ ^l)- After some algebra, we express: 

n ^ ■ (Ap + Bq) + l{Bp - Aq) 

Ok + ink = =-7^ (10) 

UqUi 



Then the Kerr angle is: 



(Ap + Bq) ^^^^ 



and the ellipticity is: 



Aq-Bp 

where A = Ra^yRa^x + la^yla^x, B = la^yRa^x - Ra^yla^^, Dq = {Ra^xY + (Icr^x)^ and 

In this representation, an explicit analytic relation for the Kerr spike can be derived using 
the Kramers-Kroning transformation: whenever the ellipticity vanishes, a spike in the Kerr 



spectra appears. [311] Thus, imposing that rjk = 0, we derive an analytic relation for 9k which 
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holds for the maximum or minimum in the Kerr spectra corresponding to zero elhpticity. 
From rjk = 0, we have either (a) A = ^ or (h) B = ^ and using Eq. [11] we obtain: 



= -aI^ = -bI±^ (13) 



A necessary condition to produce a Kerr spike is that ^ = |- 
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FIG. 1: Total and atomic- angular momentum- projected density of states (DOS) in Ga.75Mn.25 As. 
Positive and negative curves refer to spin-up and spin-down components, respectively. The Fermi 
level is fixed at zero energy. 



FIG. 2: (Color on line) Spin-polarized band structure for Ga.75Mn.25GAs along the symmetry lines 
F-X-M-F-R-M in the irreducible BZ plotted without spin-orbit interaction for clarity: parts (a) and 
(b) refer to spin- up (-down), respectively. The horizontal line marks the Fermi level (zero energy). 
Numbers in round brackets at F label groups of bands discussed in the text. The vertical lines 
indicate the different possible band to band transitions at 0.4, 1.1, 2.8, 3.8 eV within 0.1 eV (solid 
black, dotted blue, dashed red and dot-dashed orange respectively). The labels for symmetry points 
and lines are according to the standard group-theoretical notation of Bouckaert, Smoluchowski, and 
Wigner. Spin-orbit effects are neglected. 
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FIG. 3: Calculated real (a) and imaginary (b) parts of diagonal and off-diagonal components of 
the optical conductivity. The Drude term has not been included in the spectra. Vertical lines in 
Re 

Gxx indicate features discussed in the text. 



FIG. 4: Selected contributions from spin-up (solid line) and spin-down (dashed line) transitions to 
the total Re cTxx spectra. The joint density of states (JDOS) is also reported, in arbitrary units. 



FIG. 5: (a) Real part of the diagonal conductivity: our calculated spectrum (bold line) is compared 



with the results of model calculations (dashed line) from Ref. 



from Ref. 



6a] and experiments (dotted line) 



The inset in panel (a) shows the difference between the interband and the total 



(interband + intraband) conductivity. Panels (b) and (c) show the calculated real and imaginary 
parts of the dielectric constant, respectively, with (bold line) and without (dashed line) A fitting. 
Ellipsometry results from Ref. [66i] are also shown (solid line). 



FIG. 6: Prospective view of the irreducible wedge of Brillouin zone showing the origin of the 
interband transitions at ~ 3.8 eV. The length of the arrows shows the relative contributions to Re 
axx in the energy interval 3.6-4.0 eV (see Fig. [3]). Increasing of the length of the arrows represents 
an increasing contribution. See text for further details. 



FIG. 7: (Color online) (a) Ellipticity and Kerr angle; the main features of the Kerr spectra are 
indicated by arrows; (b) Different contributions to the numerator of the formula presented in 
Sect. |Vl] and the Appendix (arrows indicate the main spectral features); (c,d,e) decomposition 
of the denominator (see text for details). Arbitrary units are used except for the Kerr angle (in 
decimal degree). See text for details. 
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